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O . Abstract 

^^ I We investigate two-flavour and two-colour QCD at flnite temperature and 

^S^ • chemical potential in comparison with a corresponding Nambu and Jona- 

^^. Lasinio model. By minimizing the thermodynamic potential of the system, we 

Qh! confirm that a second order phase transition occurs at a value of the chemical 

potential equal to half the mass of the chiral Goldstone mode. For chemical 

potentials beyond this value the scalar diquarks undergo Bose condensation 

and the diquark condensate is nonzero. We evaluate the behaviour of the 

/S \ chiral condensate, the diquark condensate, the baryon charge density and the 

j^ ' masses of scalar diquark, antidiquark and pion, as functions of the chemical 

potential. Very good agreement is found with lattice QCD {Nc = 2) results. 

We also compare with a model based on leading-order chiral effective field 

theory. 
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1 Introduction 

The phase structure of QCD has been subject of intense investigations in recent 
years. Precise numerical data have become available concerning QCD thermody- 
namics at high temperature via large-scale computer simulations on the lattice (for 
a review see [1]). The study of full QCD at finite baryon density is still a formidable 
challenge, due to the limitations of standard Monte Carlo simulations when applied 
to systems at finite chemical potential (for recent results see [2, 3]). Present develop- 
ments are aimed at improved strategies [4] to deal with the fact that the determinant 
of the Euclidean Dirac operator becomes complex at finite chemical potential. 

An interesting perspective of finite-density QCD is the emergence of colour su- 
perconductivity (CSC). This was revealed first by calculations based on one-gluon 
exchange: Barrois, Bailin and Love [5, 6] and later Iwasaki and Iwado [7] pointed 
out that the induced attractive force near the Fermi surface creates quark Cooper 
pairs resulting in CSC in the case of QCD at low temperature and high density. 
In the late nineties, using an instanton model of the effective interaction, Alford, 
Rajagopal and Wilczek [8, 9] and Rapp, Schafer, Shuryak and Velkovsky [10] argued 
that the energy gap is expected to be of the order of 100 MeV. 

No first principle computations exist at this moment concerning the phenomenon 
of colour superconductivity in full Nc = 3 QCD. One response to this situation has 
been to start from simpler QCD-like theories with additional antiunitary symmetries 
that guarantee the Fermion determinant to be real at non-zero chemical potential 
and therefore allow the study of such theories on the lattice. Examples of such ex- 
plorations include QCD with two colours and fundamental quarks and QCD with an 
arbitrary number of colours and adjoint quarks [11]. The physics of both these the- 
ories is quite different from full three-colour QCD. Nevertheless these differences are 
easily understood and classified. Knowledge of the critical conditions for phase tran- 
sitions in these schematic cases may offer qualitative clues about critical phenomena 
encountered in three colour QCD, such as diquark condensation. 

In two-colour QCD, diquarks can form colour singlets which are the baryons 
of the theory. The lightest baryons and the lightest quark-antiquark excitations 
(pions) have a common mass, m.^, and this spectrum determines the properties of 
the ground state for small chemical potentials. General arguments [12] predict a 
phase transition from the vacuum to a state with finite baryon density at a critical 
chemical potential fic, which is the lowest energy per quark that can be realized 
by an excited state of the system. This state is populated by light diquarks, and 
one expects fic = itLjj/2. The Bose-Einstein condensation of diquarks, with nonzero 
baryon number, can be interpreted as baryon charge superconductivity. 

The (T, fi) phase diagram of QCD with two colours has been studied by Dagotto 
et al. using a mean- field model of the lattice action [13, 14]. The smallness of Hc 
has been exploited to study the zero temperature phase transition using a chiral 



effective Lagrangian extended to tlie flavour symmetry SU{2Nf) [11, 15-18]. Otlier 
approaclies to two-colour QCD have also been explored, based for example on a 
random matrix model [19,20] and on the renormalization group [21]. Several of 
these model calculations have been verified by lattice simulations [22-41]. 

In the present paper we investigate the relationship between Nc = 2 QCD and 
a corresponding Nambu and Jona-Lasinio (NJL) model [42-46] in which gluonic 
degrees of freedom are "integrated out" and replaced by a local four-point interaction 
of quark colour currents. This amounts to effectively replacing the local colour gauge 
symmetry by a global one, with the assumption that coloured (gluonic) excitations 
are far removed from the low-energy spectrum and hence "frozen" . Similar models 
have already been used to study the QCD colour superconductivity phase with two 
[47-51] and three flavours [52-54] (for a recent review see [55]). The specific aim 
of this work is to test the effectiveness of the NJL model, with its dynamically 
generated quasiparticles, in reproducing the thermodynamics of two-colour QCD, 
and to compare our results quantitatively with those obtained from recent lattice 
computations. We study the behaviour of the chiral and diquark condensates, and of 
the baryon density, as functions of temperature and chemical potential, both in the 
chiral limit and for finite values of the current quark masses. We investigate, again 
for both zero and finite quark masses, the two-colour QCD phase diagram in the 
T-fi plane. As further applications we evaluate the pion, diquark and antidiquark 
masses, as functions of the chemical potential. We compare our results to lattice 
data and also to the predictions from chiral effective field theory. 

2 Two colour NJL model 

Consider as a starting point the Lagrangian 

3 

C = ^{x) {tYd, - mo) ^ (x) -G,J2 "^M (^) Ja (^) , (1) 

a=l 

with a four point interaction that represents the local coupling between colour cur- 
rents J^ = '(l)'^^t""ip involving the quark fields ip and the SU{2)coiour generators ta 
with tr{t°'t^) = 2Sab- Here Gc is an effective coupling strength of dimension (length)^ 
and mo is the diagonal current quark mass matrix. 

In this paper we restrict ourselves to the case of two quark flavours {Nf = 2). 
In this case there are only two order parameters, the quark condensate (ipip) and 
the scalar diquark condensate, symbolically denoted by (ipip)- It is convenient to 
rewrite the interaction between quarks, by Fierz transformation, in terms of the 
colour singlet pseudoscalar/scalar quark-antiquark and scalar diquark channels. The 
resulting Lagrangian reads 

jC-njl = ■?/' (i7^i9^ - mo )■?/' + £gg + £gg + (colour triplet terms), (2) 



TT 

where G and H are constants which describe quark-antiquark and quark-quark 
interactions, respectively, ta are Pauh matrices in colour space and Ti are Pauli 
matrices in flavour (isospin) space. We have introduced the charge conjugation 
operator for fermions: 

C = i7o72- (3) 

The coupling constants G and H in the Lagrangian (2) are fixed by Fierz trans- 
forming the colour-current interaction in (1) to obtain 

G = H=^G, (4) 

(see the Appendix for details). 

As mentioned, the local SU{Nc = 2) gauge symmetry is replaced by global 
SU{2)coiour in this model. In the chiral limit, the Lagrangian (2) is invariant under 
an enlarged flavour symmetry SU{Nf) x SU{Nf) x U{1) -^ SU{2Nf), which con- 
nects quarks and antiquarks: the so-called Pauli-Giirsey symmetry, a characteristic 
feature of two-colour QCD. This symmetry relates pions and scalar diquarks. It is a 
natural ingredient of the "equivalent" NJL model, with eq. (4) relating the coupling 
constants of the model Lagrangian. 

Starting from the Lagrangian (2) and using standard bosonization techniques, 
we introduce the auxiliary scalar (a), pseudoscalar triplet^ (vf) and diquark (A, A*) 
fields, thus obtaining the following equivalent Lagrangian in the colour singlet sector: 

C = i) {i'^^d^ -vdo + a + ij^T ■ vf) ^ - -A>'^ 6*75X2^2^ 

+ ^^V575r2t2CV5" - ^^ - ^. (5) 

It is useful to represent the quark fields by a bispinor defined in the following way: 

Furthermore, we introduce the matrix propagator 



^Isovectors such as the pion field are denoted by 7?. 



(the inverse of the so-called Nambu-Gorkov propagator) where we have defined 

M = {mo - (r)l - i-f^T ■ n- (8) 

here 1 = Ic ■ 1/ ■ Id is the unit matrix in colour, fiavour and Dirac indices. We 
consider the fiavour-symmetric case with rriu = rrid = mo. Integrating over q{x) and 
q{x) we obtain the effective Lagrangian in terms of the auxiliary field variables a, n, 
A and A*. It can be written as: 

^eff = -^ - ^ - . / T^^Tr In (5^ (p)) . (9) 

^^ 2G 2H J (27r4)2 ^ ^^'^ ^ ' 

The trace in this expression is taken over fiavour, colour and Dirac indices, and the 
factor I compensates for double-counting in the q and q fields. 

Solving the field equations for cr, vf, A and A* and working in the mean field 
approximation^, we can evaluate their vacuum expectation values. The mean field 
value (tt) of the pseudoscalar isotriplet field is always equal to zero. The a field has 
a non-vanishing vacuum expectation value as a consequence of spontaneous chiral 
symmetry breaking, while the diquark fields A and A* are expected to have nonzero 
mean values only in dense matter. An interesting limiting situation is encountered 
when mo = (chiral limit) together with fi = 0. In this limit the extended SU{2Nf) 
symmetry with Nf = 2 (and G = H) implies that the thermodynamic potential 
depends only on i?^ = o"^ + |Ap so that there is a degeneracy along the circle with 
constant radius R. This case will be further discussed in section 5. 

After solving the field equation for a, we can work in terms of the effective quark 
mass m, which is related to (cr) through the self-consistent gap equation 

m = mo — (o-) = mo — G{ijjip). (10) 

Note that (a) = G{ipip) is negative in our representation, and {ipip) = {ipu'^u + i'd'^d) 
with {tpui^u) = {i'd'4'd)- 

3 Parameter fixing 

The three parameters of the model are the "bare" quark mass mo, a loop-momentum 
cutoff A and the coupling strength G = H. Even if we are considering the Nc = 2 
NJL model, we choose to reproduce the known chiral physics in the hadronic sector. 
This is reasonable since, in colour singlet channels, Nc enters only parametrically 
in the relevant physical constants and observables. For this reason, we fix those 
parameters through the constraints imposed by the pion decay constant and the 
chiral (quark) condensate: 

^In the mean field approximation the fields are replaced by their expectation values for which 
we will later on continue using the notation a and A for simplicity and convenience. 



A [GeV] 


G = H[GcY-^] 


mo[MeV] 


m[MeV] 


|(^„^„) 1/3 [Me V] 


/.[MeV] 


m^[MeV 


0.78 


10.3 


4.5 


361 


259 


89.6 


139.3 



Table 1: Parameter set used in this work, and the corresponding physical quantities. 



The pion decay constant /,r is evaluated in the NJL model through the follow- 
ing relation: 



.2r(i) 



(1) 



/^ = Am I^ (m) where I^ {m) = —iN^ 



(27r) (p^ — m^ + ie) 



The empirical value is f.,^ = 92.4 MeV. 
The quark condensate becomes 

{■ipu-ipu) = -Amlf^ (m) 

with 

/» (m) = iNc / T^^ ^ -. 

Its "empirical" value derived from QCD sum rules is 

(^„^„)V3 ^ (V5^^,)V3 ^ _ (240 ± 20) MeV. 



:iii 



(12) 
(13) 

(14) 



The current quark mass mo is fixed from the Gell-Mann, Oakes, Renner 
(GMOR) relation: 

2 -mo (V^V^) 



m^ 



P 

J -K 



(15) 



In the chiral limit, m^ = and m^ = 0. 



The Goldberger-Treiman relation, which determines the pion-quark coupling Qt^, 
follows from the previous relations: 



m = gnfw 



(16) 



with^2^(4/f(m))-i. 

We will first perform all our calculations with a finite value for the current quark 
mass mo, and then investigate the chiral limit, mo -^ 0. The parameters obtained 
by imposing the constraints (11-15) are shown in Table 1. 



4 Results at finite T and /i 

We now extend the NJL model to finite temperature T and chemical potentials /i 
using the Matsubara formalism. We consider the isospin symmetric case, with an 
equal number (and therefore a single chemical potential) of u and d quarks. The 
quantity to be minimized at finite temperature is the thermodynamic potential: 

n{T,^^) = -T$^y"^iTrln(^i^-^(^c.„,p-)) 

+ 2G + W' ^''^ 

where cj„ = (2n + l)7rT are the Matsubara frequencies for fermions and the inverse 
quark propagator including the chemical potential fi is now defined as 

^-1 (n^ n) = ( ^ ^ ^ ^ ^^° A75r2t2 

Using the identity 

Trln(X) = lndet(X) (19) 

we can evaluate the trace in (17) and obtain 

^-'..„„^V4,„fifi±r)!).41„(ii±r-l. ,20, 




where we have defined E^ = -\/(e^) + |A| , with e^ = e ± /i, e = y^p ^ + m^. Next 
we evaluate the Matsubara sum in eq. (17) using the following relation: 

^ E 1^ y2 = i^^ + 2Tln (1 + exp {-E^/T)) . (21) 

n=—oo \ / 

The thermodynamic potential becomes: 
niT,^.) = -4|^[2Tln(l + exp(-^))+ (22) 



+ 2TlnM +exp(-:^j j + (E+ + E-) 



2G 2H 



In eqs. (20)- (22), the effective (constituent) quark mass m is related to the current 
quark mass and the a field through eq. (10). 

The mean values for the a and A fields are determined by minimizing the ther- 
modynamic potential. One obtains the following set of coupled equations that must 
be solved simultaneously in order to find the solutions for a and |A|: 
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Figure 1: Scaled expectation values {a) and (|A|) as a function of the chemical potential for 
different temperatures. Continuous lines correspond to T = 0, dashed lines to T = fOO MeV, 
dotted lines to T = 150 MeV and the dashed-dotted line corresponds to T = 200 MeV ((|A|) =0 
in this case). 



a 



^2 / dPP 



e — n e + n 



E^ 



(23) 



e — ji 



e-\- jji 



(exp (^) + 1) E~ (exp (^) + l) E- 



lAI 



-T I ^PP^ 



TX^ 






lAI 



lAI 



lAI 



E^ Uexp(^) + 1)E- (exp (^) + 1) i?+ 



In Fig. 1 we show our results for the scaled expectation values of the a and A fields 
as a function of the chemical potential for different temperatures. One observes that 
at T = the system undergoes a second order phase transition at a critical chemical 
potential [ic = fn^/^, as predicted by general arguments. The value of the pion 
mass that we consider here is the one evaluated in the model and shown in Table 1. 
So this model exhibits diquark condensation at chemical potentials larger than fic, 
where the value of the chiral condensate is correspondingly reduced. At T = 0, 
A is always non-vanishing for /i > fic'- the diquark phase persists for large fi. For 
temperatures T > 200 MeV, on the other hand, the diquark condensate vanishes 
even for large chemical potentials. 

The chiral effective Lagrangian approach [11] predicts the following behaviour 



for the diquark condensate as a function of the chemical potential at /i > fie'- 



{^ (|A|) /. (m^ , ^24) 



l(V'V^)o| IWol V V2/i 

which means that (|A|) should reach the vacuum expectation value of the (scaled) 
chiral condensate asymptotically as /i — *> oo. In the NJL model, the scale of variation 
for ji > He IS set by the momentum cutoff A. As a consequence, |A(/i)| increases 
until /i ~ A (corresponding to /i/m^ ~ 5). For larger values of ^ the relevant 
interactions become weaker and |A| tends to decrease with /i. This feature is an 
artifact, however, since the applicability of the NJL model is limited to energy 
and momentum scales below A. For chemical potentials smaller than the cutoff 
scale the agreement between NJL and chiral Lagrangian calculations is excellent, as 
expected. At very large chemical potential, perturbative gluon exchange presumably 
takes over, with decreasing interaction strength as /i increases. 

In Fig. 2 we show a comparison of our results for the scaled chiral and diquark 
condensates at T = as a function of the chemical potential, with lattice data taken 
from ref. [27]. These data have been obtained by studying two-colour QCD with 
staggered fermions in the adjoint representation. It was found that the positive 
determinant sector behaves like a two-flavour theory. As we can see, the agreement 
of our results with lattice data is remarkable. The dashed lines are the predictions 
from chiral effective field theory. 

In Fig. 3 we show the scaled (a) and (|A|) as a function of the temperature for 
different values of the chemical potential. In this way we find, as a function of the 
chemical potential, the critical temperature of the phase transition, so that we can 
draw the phase diagram of two-colour QCD as modelled in the NJL model. We 
show it in Fig. 4. At very small chemical potentials we have a transition from a 
system in which chiral symmetry is spontaneously broken to a system where it is 
restored (from region I to region II) with (|A|) = in both phases. Region III is 
the superfluid phase with (|A|) ^ 0. 

An interesting quantity is the baryonic density 

p^-^. ,25) 

The lattice data of ref. [27] show a scaled baryonic density defined as: 

Leading order chiral effective field theory [11] gives the following behaviour at /i > ^c- 

^ ' ^ ' ■'■^" ' ' (27) 



2m,r 





0.25 0.5 0.75 1 1.25 1.5 1.75 
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Figure 2: Scaled (a) and (|A|) as a function of the chemical potential at T = 0: our results 
(solid lines) are compared to the lattice data taken from ref. [27]. The different symbols (open 
circles, squares and diamonds) for the chiral condensate correspond to different values for the quark 
masses. The dashed lines are the predictions from chiral effective field theory [11]. 
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Figure 3: Scaled (a) (a) and (|A|) (b) as a function of temperature for different values of [ijmT^. 
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Figure 4: Phase diagram in the NJL model with two colours. The zone I is a region in which 
chiral symmetry is spontaneously broken, and (|A|) = 0; in region II chiral symmetry is restored, 
and again (|A|) — 0; region III is the superfluid phase in which (|A|) ^ 0. 

Fig. 5 presents our results for the scaled baryonic density (26) as a function of 
the chemical potential at zero temperature, in comparison with the lattice data 
for the same quantity. Our results are in good agreement with lattice data at 
moderate chemical potentials, while for large chemical potentials the baryon density 
is underestimated. This difference may be caused by the mean-field approximation. 
Correlations between quasiparticles, not covered by this approximation, tend to 
become increasingly important with growing density. 

4.1 Pion and scalar diquark properties 

This Section presents our results for the masses of the (pseudo)-Goldstone bosons, 
namely the pion, the scalar diquark and the corresponding antidiquark. 

In order to evaluate the masses of the bosonic fields, we expand the effective 

action 

AA* 



S, 



eff 



d X 



2G 



2H 



I 

-tr 

2 



d^a;ln(5'i(a;)). 



(2J 



in a power series of the meson and diquark fields around their mean field values. 
The second-order term of this expansion identifies the mass spectrum of mesons and 
diquarks. The resulting effective action in momentum space has the following form: 



5g)^(a,7f,A,A*) 



a 



TX 



2G 



AA* i 

h -tr 

2H A 






SoASoA 



(29) 
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Figure 5: Scaled baryonic density as a function of the chemical potential at T = (Continuous 
line). The lattice data are taken from rcf. [27]. The different symbols correspond to different values 
for the quark masses. The dashed line is the prediction from chiral effective field theory [11]. 

where 5*0 is the Nambu-Gorkov propagator (18) evaluated at the mean field values 
for the bosonic fields, and A is a matrix defined in the following way: 



A 



a + 2757? ■ f 
-A*75r2t2 



A75r2t2 
a — i'js^n ■ r 



(30) 



(see also [56]). By analyzing the second order action (29), one observes that mixing 
terms arise, at /i > fic, between the a, A and A* fields: these terms are proportional 
to |A|, and the mixing occurs because the presence of a nonzero diquark condensate 
spontaneously breaks the baryon number symmetry. This feature was already found 
in [39]. The mass matrix turns out to have the following form: 



/ 925, 



M 



(2) 
eff 







d^€\ 







dcrdA dadA* 



\ 



\ 



dAda 
dA'du 9A*9A 



9A2 



0^ 

aA9A* 
aA*2 



(31) 



/ 



and the masses of the various modes are found by solving the equation 



det (M) = 0. 



(32) 
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Figure 6: Scaled pion mass as a function of ii/m^r at T = (continuous line). The lattice data 
are taken from ref. [27] and have been rescaled in order to show dimensionless quantities. The 
different symbols correspond to different values for the quark masses. The dashed line is itIt^ — 2/i, 
as predicted in leading-order chiral effective field theory [11]. Also shown is the (scaled) pion decay 
constant /vr/TiTr and its evolution with increasing /i. 

Evidently the pion fields do not mix with the others, while the a, the diquark and 
the antidiquark fields mix in the phase with |A| ^ 0. 

The behaviour of the scaled pion mass as a function of the chemical potential 
is shown in Fig. 6, in comparison to the lattice data. The pion mass increases 
linearily with the chemical potential at /i > /ic- This behaviour was anticipated in 
the calculations by Kogut et al. [11]. They in fact predicted for rriT, the following 
behaviour at /i > ^c- 

rriT, = 2/i, (33) 

as indicated by the dashed line in Fig. 6. Our result is in very good agreement with 
both the lattice data and the predictions using the leading-order chiral effective 
Lagrangian. 

The behaviour of the pion and a masses and of the pion decay constant as 
functions of temperature at /i = is shown in Fig. 7. At temperatures T exceeding 
the critical Tc for the chiral transition at which (a) tends to zero, m„ becomes equal 
to the pion mass and both masses rise continuously with increasing T. The pion 
decay constant tends to zero at the same time. 

Next, consider the other two bosonic modes of the theory: the scalar diquark 
and its antidiquark. The behaviour of their masses at finite chemical potential is 
shown in Fig. 8 in comparison to the pion mass: at /i = they are all degenerate, 
as predicted on the basis of general arguments, but they behave in different ways as 
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Figure 7: Pion mass, (7 boson mass and pion decay constant as a function of temperature at 
/x = 0. 



the chemical potential increases: for fi < fi^ 



m 



(0) 



/2 the pion, which does not carry 



baryon charge, is not affected by fi, while the diquark and antidiquark masses are 
shifted according to their baryon number B = ±1. They follow in fact the behaviour 
observed also in chiral effective field theory [f f]: 



rriA = itLt, — 2fi, 



rriA* = rriT, + 2/i. 



(34) 



For fi > fic, the appearance of the diquark condensate spontaneously breaks the 
baryon number symmetry. The scalar modes (diquark, antidiquark and sigma) get 
mixed. The new eigenmodes are linear combinations of the original quasiparticle 
states. By solving eq. (32) we find the masses of the new orthogonal modes. One 
of them, which we denote by A, is massless and can be identified with the true 
Goldstone boson of the theory, corresponding to the spontaneous breaking of the 
baryon number (f/(l)) symmetry. The other two modes are massive. One of them, 
which we denote by A*, follows the behaviour derived in the paper by Kogut et al: 



Ma. 



2/i-y/l + 3(m„/2/i)' 



(35) 



5 Chiral limit 

In the chiral limit mg — > (m^ -^ 0), and at ;U = 0, the thermodynamic poten- 
tial (22) (with G = H) is a function only of cr^-|-|Ap, as already mentioned. This is a 
natural outcome once the relation between the coefficients G and H is fixed through 
the Fierz transformation of the colour current-current interaction (see eq. (4)). As 
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Figure 8: Spectrum of pions and diquarks/antidiquarks as a function of the (scaled) chemical 
potential at zero temperature. 

a result, Q is invariant under the rotation which connects the chiral and the diquark 
condensate along the circle a^ + |Ap=const. Because of this symmetry, the chiral 
condensate is indistinguishable from the diquark condensate for mo = /x = 0, so 
that a state with finite (a) can always be transformed into a state with finite (|A|) 
and (a) = 0. The phases with spontaneously broken chiral and baryon number 
symmetries are degenerate in this limit. 

As soon as the chemical potential takes a finite value, the favourable phase is 
the one with a non-zero diquark condensate and zero chiral condensate. This is 
evident from Fig. 9 which shows the contour plots of the thermodynamic potential 
as a function of a and |A|. In the left panel we have T = fi = and the rotational 
invariance is evident. In the right panel we have introduced a very small chemical 
potential, which is nevertheless sufficient to break the rotational invariance along 
-R^ = (a)^ + (|A|)^ and favour the phase in which (a) = and (|A|) ^ 0. 

Minimizing the thermodynamic potential of the system, one finds the mean-field 
values of the chiral and diquark condensates. Our results in Fig. 10 display (|A|) as 
a function of temperature for different chemical potentials. The chiral condensate 
is always equal to zero in those cases. 

In Fig. 11 the phase diagram of the two-colour NJL model in the chiral limit is 
compared to the one using a finite value of the bare quark mass rriQ. As one can see, 
the phase boundaries for rriQ = and rriQ ^ become identical at large chemical 
potentials, whereas at small /i they show a qualitatively different behaviour. In the 
exact chiral limit there are only two phases in the theory: the superfiuid phase with 
(|A|) 7^ and the high-temperature phase with (|A|) = 0, separated by a critical 
temperature of about 0.2 GeV. 
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Figure 9: Contour plots of the thermodynamic potential in the chiral limit (too = 0) as a function 
of cr and A for T = /x = (a) and T = and ^i = 20 MeV (b). 



Consider next the pion and diquark masses in the chiral hmit and their variations 
with increasing chemical potential. The chiral condensate is always equal to zero 
in this limit. Consequently, the A mode is a true Goldstone boson and its mass 
is always equal to zero, while the A* and pion masses are degenerate. Explicit 
symmetry breaking by a finite chemical potential lets these masses scale as m^, = 
rrij^ = 2/i. The degeneracy of A* and tt is removed as soon as a small non-zero quark 
mass rriQ is introduced. This also gives a finite mass to the A mode, which is again 
equal to zero above fj,c = m^(o)/2. 

Fig. 12 illustrates this situation for a very small value of itlq (~ 0.1 MeV). The 



critical value fic of the chemical potential is identified as fie 



m 



(0) 



/2, as discussed 



previously, but now of course with a very small value of the vacuum pion mass rrin . 
As the limit uiq ^ is approached, mi -^ and /ic ^ 0: the low-temperature 
system is always in the superfiuid phase for any value of /i. At /i = we recover the 
exact Pauli-Giirsey symmetry, with vanishing pion and diquark masses. 



6 Conclusions 

We have investigated a two-colour and two-fiavour Nambu and Jona-Lasinio model 
at finite temperature and finite baryon chemical potential, with the primary aim 
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Figure 10: Mean field value of the |A| field as a function of temperature for /i = (continuous) 
and ^ = 350 MeV (dashed). 
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Figure 11: Comparison between the phase diagram of two-colour QCD in the chiral limit (con- 
tinuous) and for bare quark mass mo ^ (dashed). 
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In the exact chiral hmit, m 



A* 



2/i and m^ = 0. 



of exploring the capability of such a model to reproduce the thermodynamics of 
Nc = 2 lattice QCD. The starting point is the assumption that gluon dynamics can 
be integrated out and reduced to a local interaction between quark colour currents. 
By Fierz rearrangement, this implies a one-to-one correspondence between inter- 
actions in colour singlet quark-antiquark and diquark channels (the Pauli-Giirsey 
symmetry). 

The resulting spontaneous (dynamical) symmetry breaking pattern identifies 
pseudoscalar Goldstone bosons (pions) and scalar diquarks as the thermodynam- 
ically active quasiparticles. The successful comparison with Nc = 2 lattice data 
indicates that this simple NJL model does indeed draw a remarkably realistic pic- 
ture of the quasiparticle dynamics emerging from Nc = 2 QCD, even though the 
original local colour gauge symmetry of QCD has been reduced to a global colour 
SU{2) symmetry in the NJL quasi-particle model. We note that colour (triplet) 
quark-antiquark modes which are the remnants of gluon degrees of freedom in this 
model, are far removed from the low-energy spectrum. Poles of the respective Bethe- 
Salpeter amplitudes appear at mass scales several times the NJL cutoff scale [57]. 

We confirm that a diquark condensate develops at chemical potentials /i > yU^ = 
m7r/2. The correlated evolution of the chiral and diquark condensates with increas- 
ing /i, as observed in A'"c = 2 lattice QCD, is very well reproduced. Had we started 
from NJL four-point interactions with independent, arbitrary coupling strengths 
in quark-antiquark and diquark channels, the condensate pattern would have been 
quite different. It appears that modelling the low-energy dynamics of A^c = 2 QCD 
is already done surprisingly well when using just a colour current-current interaction 
with a single strength parameter. 

The calculated baryon density, obtained by taking the derivative of the thermo- 
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dynamic potential with respect to the chemical potential, describes the correspond- 
ing lattice results well in the range /i < 2/1^. Deviations occur at larger /i which 
presumably indicate the increasing importance of correlations between quasiparticles 
beyond the mean-field approximation. 

The NJL model also permits an instructive study of the way in which this system 
behaves in the chiral limit which is not directly accessible in lattice computations. 
In particular, the limits of vanishing quark mass and vanishing baryon chemical 
potential do not commute, as expected, and have to be handled with care. 

The low-energy physics of QCD differs qualitatively between Nc = 2 and Nc = 3 
because of the very different nature of the baryonic quasiparticles in these two the- 
ories. Nevertheless, the success of the present studies encourages further extended 
investigations also for Nc = 3 thermodynamics, using NJL type quasiparticle ap- 
proaches above the critical temperature for deconfinement, in close contact with 
lattice QCD simulations. 



We thank Pietro Faccioli, Jiri Hosek, Georges Ripka and Michael Thaler for stimu- 
lating discussions and valuable comments. 

Appendix 

We start from the colour current interaction (1) and show that performing a Fierz 
transformation we obtain the Lagrangian (2) with the coupling coefficients related 
by (4). 

In order to demonstrate this identity for the coefficient of the scalar diquark 
interaction (H) we must Fierz transform this interaction into the qq channel, while 
for G we must Fierz transform into the qq channel. 

Let us start with H; we rewrite the interaction term of eq. (1) and keep track 
explicitly of all colour, flavour and Dirac indices: 



^L = -aX](V^7«t»' (36) 

a=l 
3 

= -Gc Y^ m,p,p,i^j,q,u^k,r,p^l,s,a (7a)^^ (t")^^ {Qij {Q kl ^Pq^r 



with: 



a=l 



i,j,k,l — > colour indices 

p,q,r,s — > flavour indices 

/i, I', p,(T — > Dirac indices. 
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We start by performing the Fierz transformation for the flavour indices using the 
following relation 



1 

^pq^rs = T;"^ in) pr in) sg, 



6=0 



where we have defined 

1 
1 



To 



and Tb = Pauli matrices with b = 1, 2, 3 



)■")'-' ) 



(37) 



(3J 



thus obtaining 



1 



3 3 



^L == -2'^cJ2Yl [^i,P,l^'^J,1,'^^k,r,pAs,cr ila)f,^ (t")^. (Qij (Qkl in)pr (u) ^^^ (39) 
a=l 6=0 

In order to Fierz-transform the colour indices we use the relation 

3 



J2^^-^^A^'^' 



1 



kl 



a=l 



S^kSlJ + ih\k itl)l, + (^3),, {h\^\ - - (t2),fc {h\^ , (40) 

thus obtaining 

1 ^ 



■'int 



6=0 



+ (^i)ifc (^i)«i + ih)ik {h)ij - 3 (t2)ifc (t2)jjj (r6)p, (rft)^^ 
At the end we perform the Fierz transformation for the Dirac indices and find 

^ (7.C*)^, (C7"),. + 



(41) 



'-'mt 



1 ^ r 



6=0 



^i,p,M^i,g,^^fc,r,p^/,s,a ( (C*)^^ (C)^^ 2 



- \ bal^C*)^^ {Cri5)au + in,C*),, (^C75).. j {S^kS^j + 

+ itl)^k (^l)/j + (^3)ifc (^3)zj - 3 (t2),fc (t2)zj) (^6)^, (u) ^^ = 
1 3 _ _ 



4 



6=0 5=0,1,3 



1 - - 1 - - 

6=0 



(42) 



20 



where we have introduced the charge conjugation matrix operator for fermions C = 
27072- We can easily read from eq. (42) the coefficient of the scalar diquark channel: 



H = -Gr. 



(43) 



Next we show that also G = 3Gc/2, starting from eq. (36) and performing a Fierz 
transformation into the qq channel. 

We start from the flavour-S'[/(2) identity 

(44) 



SpgSrs = 2 5Z (^^) 
6=0 



ps \ "frq 



and obtain 



3 3 



'-'mt 



2^^ Y^ Y^ [Ap,f^^j,g,'^^k,r,pi^i,s,- (7«)/,^ (7")p^ {ta)ij {ta)ki {rb)^^ in)^^ .(45) 



a=l b=0 

Then we transform colour indices by using 

3 



Y^ {ta)ij {ta)kl = T^^ilhj -T^Yl (^^)i« (^^ 



'fej 



(46) 



a=l 



c=l 



and find 



C 



int 



1 ^ 

jGc ^ m,p,^,i'j,q,u^Kr,pi^Us,a (7a)^^ (7°)p<, (3(5i,4j 



6=0 



Y (^-)*^ (^-)fcj (^^)p« (^^)rg 



c=l 



(47) 



Finally the Dirac Fierz transformation leads to: 

3 

I 
C 



■'int 



lG=E 



6=0 



V'i,p,M^j,g,^^fc,r,p^/,s,a ( VV - 2 '^^"''^'^ '^^"'''^^ ~ 2 '^^"^^''a''^ *^^"^^'' 



pi^ 



+ (^75)^^ (^75)pJ 35ii4i - XI (^^)i« (^^ 



'fcj I V''"bJps Vbjrq 



c=l 



i«=i:|3 



6=0 



(^Tb^) + {ii)'^r>ni)) - - (ip'jarb^) ~-{i'iai5n^y 



3 r 



- E 



o=l 



{iptaTblp) + [ii^l^taUi)) - -{ipl^taUi)) - -{i^l,j,l-,tani))' 



from which we can easily read 



Gr — — Gr- 

2 



(4J 
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